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Abstract 

We present the initial results of multi-epoch VLBI observations of the 22 GHz H2O masers in the 
Orion KL region with VERA (VLBI Exploration of Radio Astrometry). With the VERA dual-beam 
receiving system, we have carried out phase-referencing VLBI astrometry and successfully detected an 
annual parallax of Orion KL to be 2.29±0.10 mas, corresponding to the distance of 437±19 pc from the 
Sun. The distance to Orion KL is determined for the first time with the annual parallax method in these 
observations. Although this value is consistent with that of the previously reported, 480±80 pc, which is 
estimated from the statistical parallax method using proper motions and radial velocities of the H2O maser 
features, our new results provide the much more accurate value with an uncertainty of only 4%. In addition 
to the annual parallax, we have detected an absolute proper motion of the maser feature, suggesting an 
outflow motion powered by the radio source I along with the systematic motion of source I itself. 

Key words: Astrometry: — ISM: individual (Orion KL) — masers (II2O) — radio lines: ISM — ISM: 
jets and outflows 



1. Introduction 

Distance is one of the most fundamental parameters in 
astronomy. However, it has been difficult to measure accu- 
rate distances to stars, galaxies, and other astronomical 
objects without assumptions. The most reliable way to 



determine the distance is an annual trigonometric par- 
allax method, based on precise measurements of posi- 
tion and motion of the object. In 1990's, the Hipparcos 
satellite extensively measured annual parallaxes for more 
than 100 000 stars with a typical precision of 1 mas 
level (Ferryman et al. 1995, 1997), which allowed us 
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to refine various fields of astronomy and astrophysics. 
Nevertheless, the distances measured with Hipparcos were 
limited only within a few hundred pc from the Sun, which 
was far smaller than the size of the Galaxy, 15 kpc in 
radius. 

In the last decade, phase-referencing VLBI astrometry 
has been developed, with which the position of a tar- 
get source is measured with respect to a reference source 
(Beasley & Conway 1995). Using extragalactic radio 
sources as the position references (e.g. sources listed in 
the ICRF catalog; Ma et al. 1998), we can measure the 
absolute position of the target source, which lead us to 
derive its annual parallax. With recent highly precise 
VLBI astrometry, annual parallaxes have been success- 
fully measured for the Galactic CII3OII maser sources at 
the 12 GHz band (Xu et al. 2006) and H2O maser sources 
at the 22 GHz band (Kurayama et al. 2005; Hachisuka 
et al. 2006) with the NRAO Very Long Basehne Array 
(VLB A). The annual parallax measurements with VLBI 
have also been carried out for non-thermal radio contin- 
uum emission from young stellar objects (e.g. Lcstradc 
et al. 1999; Loinard et al. 2005). The highest accu- 
racy of these VLBI astrometry is reported to be 0.05 mas, 
which provides a powerful tool to measure annual par- 
allaxes with the accuracy by two orders of magnitude 
higher than that of the Hipparcos satellite, allowing us 
to measure the distances of maser sources up to 2 kpc 
away from the Sun (Kurayama et al. 2005; Xu et al. 
2006; Hachisuka et al. 2006). 

In order to extend the VLBI astrometry of maser 
sources to the whole region of the Galaxy, we have con- 
structed a new VLBI network in Japan called VERA, 
VLBI Exploration of Radio Astrometry (Kobayashi et 
al. 2003), which is the first VLBI array dedicated 
to phase-referencing observations. Each VERA antenna 
is equipped with a unique dual beam receiving system 
(Kawaguchi et al. 2000; Honma et al. 2003), which 
enables us to observe the target and reference sources 
within 2.2 degrees separation on the sky simultaneously, 
thus facilitating more efficient phase-referencing VLBI ob- 
servations compared with the conventional fast-switching 
observations. Very recently, the first results of astrom- 
etry with VERA have been reported (e.g. Honma et 
al. 2007; Sato et al. 2007), demonstrating its high 
capability of annual parallax and absolute proper motion 
measurements. The main goal of the VERA project is 
to reveal 3-dimensional Galactic structure and kinematics 
based on the accurate astrometry of hundreds of H2O (at 
the 22 GHz band) and SiO (at the 43 GHz band) maser 
sources in the Galactic star-forming regions and late-type 
stars with the highest accuracy of 10 iias level (Kobayashi 
et al. 2003; Honma et al. 2000). 

In this paper, we present the initial results of the annual 
parallax measurements of Orion KL. Because Orion KL is 
the nearest high-mass star-forming region located at an 
estimated distance of only 480 pc from the Sun (Gcnzel 
et al. 1981), it has been recognized as one of the most 
important objects to study high-mass star-formation pro- 
cesses (e.g. Genzel & Stutzki 1989). Along with its prox- 



imity to the Sun, Orion KL is known to be one of the 
brightest H2O maser sources in the Galaxy, and hence, 
it is the best test bench for the first stage of the annual 
parallax measurements with VERA. 

2. Observations and Data Analyses 

Observations of H2O masers (6i6-523 , 22 2 3 5.080 MHz) 
in Orion KL were carried out in 19 observing sessions from 
Jan. 2004 to Jul. 2006 with VERA. In this paper, we cm- 
ployed the results of total 16 observing sessions which were 
carried out under relatively good weather conditions. A 
typical interval of observations was 1 month, while some of 
them, especially in the summer season, were a few months. 
All the 4 stations of VERA were used in most of the ob- 
serving sessions, while only 3 stations were used in part of 
the sessions (2004/027, 2004/272, and 2004/333; hereafter 
an observing session is denoted by year/day of the year). 
The maximum baseline length was 2270 km (see Fig.l of 
Petrov et al. 2007) and the typical synthesized beam size 
(FWMH) was 1.5 mas x 0.8 mas with a position angle of 
—30 degrees. 

All the observations were made in the dual 
beam mode; Orion KL and an ICRF source 
J0541-0541 (a(J2000) =05h41m38.083385s, 

S{J2000) = -05d41'49.42839"; Ma et al. 1998; 
Petrov et al. 2007) were observed simultaneously. 
The separation angle between them was 1.62 degrees. 
J0541— 0541 was detected fringes with a flux density of 
about 500 mjy in all the observations, which was suitable 
as a phase reference source. The instrumental phase 
difference between the two beams was measured in real 
time during the observations, using the correlated data 
of the random signal from artificial noise sources injected 
into two beams at each station (Kawaguchi et al. 2000). 
The typical value of the phase drift between the two 
beams was 3 degrees per hour. These results were used 
for calibrating instrumental effects in the observed phase 
difference between the two sources. 

Left-handed circular polarization was received and sam- 
pled with 2-bit quantization, and filtered using the VERA 
digital filter unit (Iguchi et al. 2005). The data were 
recorded onto magnetic tapes at a rate of 1024 Mbps, 
providing a total bandwidth of 256 MHz in which one IF 
channel and the rest of 15 IF channels with 16 MHz band- 
width each were assigned to Orion KL and J0541— 0541, 
respectively. In the earlier eight observing sessions from 
2004/203 to 2005/144, wc used the recording system at a 
rate of 128 Mbps, with two IF channels of 16 MHz band- 
width each for both Orion KL and J0541-0541. A bright 
continuum source, J0530-I-1331, was observed every 1-2 
hours for bandpass and delay calibration. System temper- 
atures including atmospheric attenuation were measured 
with the chopper-wheel method (Ulich & Haas 1976) to 
be 100-600 K, depending on weather conditions and eleva- 
tion angle of the observed sources. The aperture efficien- 
cies of the antennas ranged from 45 to 52% depending on 
the stations. A variation of the aperture efficiency of each 
antenna as a function of elevation angle was confirmed 
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to be less than 10% even at the lowest elevation in the 
observations (^^20 degrees). 

Correlation processing was carried out on the Mitaka 
FX correlator (Chikada et al. 1991) located at the NAOJ 
Mitaka campus. For H2O maser lines, a spectral reso- 
lution was set to be 15.625 kHz, corresponding to the 
velocity resolution of 0.21 km s~^. The effective veloc- 
ity coverage for the H2O maser lines, which was common 
for all the observing sessions, was ±40 km s~^ relative 
to the systemic velocity of Orion KL, an LSR velocity of 
8 km s-i. 

Calibration and imaging were performed using the 
NRAO Astronomical Image Processing System (AIPS). 
At first, amplitude and bandpass calibration were 
done for each target (Orion KL) and reference source 
(JG541— 0541) independently. Then fringe fitting was 
made with the AIPS task FRING on the phase reference 
source (J0541— 0541), and the phase solutions were ap- 
plied to the target source (Orion KL). In addition, we 
adopted the results of dual-beam phase calibration mea- 
surements as described above (Kawaguchi et al. 2000). 
Because the a priori delay model applied in the correlation 
processing was not accurate enough for precise astrome- 
try, we calibrated the visibility phase using the more ac- 
curate delay model, based on the recent achievements of 
geodynamics (Honma et al. 2007) in the analyses. In 
this model, we calibrated the fluctuation of the visibil- 
ity phase caused by the Earth's atmosphere based on the 
GPS measurements of the atmospheric zenith delay due 
to the tropospheric water vapor. 

The synthesized images were made using the AIPS task 
IMAGR with natural weighting. Even after the phase 
calibrations described above, we found that the dynamic 
range of the phase-referenced images was not high enough, 
possibly due to a residual in the atmospheric zenith delay, 
as pointed out by Honma et al. (2007). To improve 
the quality of these images, we estimated the atmospheric 
zenith delay residual as a constant offset for each station, 
which maximized the coherence of the resultant phase- 
referenced image. The atmospheric zenith delay residual 
was derived to be 0-10 cm on average, depending on the 
weather conditions, while it exceeded 20 cm in the worst 
case. As a result of this calibration, the dynamic range of 
each phase-referenced image was increased by a factor of 
up to 1.5. 

3. Results 

Figure 1 shows the cross power spectra of the H2O 
masers toward Orion KL. The H2O maser lines were de- 
tected within the LSR velocity range from —10 km s~^ to 
40 km s~^. We could not find high- velocity components 
in the LSR velocity of > 40 km s~^ and < —10 km s^^ 
(Genzel et al. 1981) possibly due to our narrower effec- 
tive velocity coverage (from —32 to 48 km s~^) and lower 
sensitivity. 

In order to reveal the overall distribution of the H2O 
masers, we at first mapped the H2O maser features in the 
Orion KL region at one of the observed sessions, 2005/081, 
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Fig. 1. Examples of scalar-averaged cross power spectra of 
Orion KL observed with the VERA Mizusawa-Iriki baseline 
(1267 km). A dashed line indicates the central velocity of the 
maser feature adopted for the parallax measurement in this 
paper at the LSR velocity of 25 km s~^. 

by the method adopted in usual single-beam VLBI obser- 
vations. The H2O maser features arc found to be extended 
over the 20" x30" region as shown in Figure 2. The distri- 
bution of H2O maser features is in good agreement with 
those in Genzel et al. (1981) and Gaume et al. (1998). 
The number of H2O maser features near source I, which 
is proposed to be a powering source of the outflow and 
the H2O masers (Menten & Reid 1995; Grcenhill et al. 
1998), is smaller than that of the results of the NRAO 
Very Large Array (VLA) observations reported by Gaume 
et al. (1998). This is because most of the maser features 
near source I are resolved out with the synthesized beam 
of VERA, implying that their sizes are larger than a few 
mas (Genzel et al. 1981; Gaume et al. 1998). 

Based on the H2O maser map at the epoch of 2005/081, 
we searched for intense H2O maser features whose cross 
power spectra observed with the Mizusawa-Iriki baseline 
(1267 km; see Fig.l of Petrov et al. 2007) were detected 
with a signal to noise ratio larger than 10 at all the 16 
observing epochs. We found that 10 maser features sat- 
isfied this criterion. Among them, we analyzed the data 
for one of the maser features at the LSR velocity of about 
25 km s~^ , which was redshifted relative to that of the sys- 
temic velocity of Orion KL, an LSR velocity of 8 km s~^, 
showing relatively less significant spatial structure in the 
synthesized images and the closure phases during all the 
observing sessions. Since the peak velocity of the maser 
feature was shifted systematically from 25.7 km s~^ to 
24.5 km s~^ during the observing period of 2 years, we 
made images of maser spots for all the spectral channels 
within the velocity range of 24.5-25.7 km s~^. and deter- 
mined the position of the maser feature taking that of the 
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Table 1. Results of the least-squares analysis for the 
annual parallax and proper motion measurements 



Parameter 



Best fit value 



5 0-5 
RA [arcsec] 



10 



Fig. 2. Distribution of H2O masers in Orion KL observed 
in the epoch of 2005/081. Crosses represent the positions of 
individual or groups of H2O maser features. Filled squares 
and circle indicate the positions of radio continuum sources 
(Gomez et al. 2005) and the maser feature analyzed in this 
paper at the LSR velocity of 24.5-25.7 km s~^. Bold arrows 
indicate the absolute proper motion vectors based on our 
study and Rodriguez et al. (2005), while a dashed arrow 
shows the proper motion of the maser feature with respect to 
source I (see text). The position offsets are with respect to 
the reference position (o(J2000) =05h35ml4.363600s, 
5(J2000) = -05d22'38.30100"). 



peak velocity channel. Although we cannot rule out the 
possibility of acceleration of this maser feature, the ob- 
served velocity shift would imply a variation of its source 
structure. Along with the velocity shift, the flux density 
of the maser feature was also highly variable as shown in 
Figure 1. The variation of the maser feature suggested by 
the velocity shift and the flux variability would affect the 
accuracy of astrometry even if the maser is bright, rela- 
tively stable, and with less significant spatial structure, as 
described later. Detailed analyses for all the H2O maser 
features will be reported in a forthcoming paper. 

Figure 3 gives the results of our position measurements 
of the H2O maser feature. As shown in Figure 3, we have 
successfully measured the movement of the H2O maser 
feature for longer than 2 years. The movement signifi- 
cantly deviates from a linear motion showing a sinusoidal 
modulation with a period of 1 year. This is clearly due 
to the annual parallax of the maser feature. In fact, the 
date of each peak in the sinusoidal curve is almost con- 
sistent with those predicted from the annual parallax of 



TT 
Ma 



2.29(0.10) mas 
2.77(0.09) mas yr^^ 
-8.97(0.21) mas yr"! 

0.36 mas 

0.74 mas 



Note — Numbers in parenthesis represent the 
estimated uncertainties. Annual parallax tt is 
derived from the right ascension data only. 

Orion KL. Assuming that the movement of the maser fea- 
ture is the sum of linear motion and the annual parallax, 
we can obtain the proper motion in right ascension /Zq 
and declination ^5, the initial position in right ascension 
ao and declination 5q, and the annual parallax tt for the 
maser feature by a least-squares analysis. 

Initially, we determined these 5 parameters simultane- 
ously, using both right ascension and declination data. In 
this case, the derived annual parallax was 2.25±0.21 mas, 
corresponding to the distance of 445±42 pc, and the stan- 
dard deviations of the least-squares analysis in right ascen- 
sion a a and in declination as were 0.36 mas and 0.74 mas, 
respectively. The larger standard deviation in declina- 
tion suggests that the astrometric accuracy in the decli- 
nation is significantly worse than that in the right ascen- 
sion. This trend can be seen in other observations with 
VERA (Honma et al. 2007; Sato et al. 2007). One of the 
possible reasons for this is that the residual of the atmo- 
spheric zenith delay would affect the astrometric accuracy, 
as discussed later. Therefore, we at first determined the 
absolute proper motion and initial position in right ascen- 
sion together with the annual parallax using the data for 
right ascension only. As a result, we obtained the annual 
parallax with higher precision to be 2.29±0.10 mas, corre- 
sponding to the distance of 437±19 pc. After the annual 
parallax was derived from the right ascension data, we es- 
timated the absolute proper motion and initial position in 
declination using the data for declination. The results are 
summarized in Table 1. 

4. Discussions 

4-.1. Astrometric error sources 

In this paper, we successfully measured the annual par- 
allax of Orion KL to be 2.29±0.10 mas through the 2- 
year monitoring observations of the H2O maser feature 
with VERA. The sinusoidal curve of the movement of the 
maser feature as shown in Figure 3 is almost coincident 
with the predicted annual parallax of Orion KL both in 
period (1 year) and phase (date of the peaks in the sinu- 
soidal curve). Therefore, the deviation from the best fit 
model, which is the combination of annual parallax and 
linear proper motion of the maser feature, should be re- 
garded as astrometric errors in our observations, rather 
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Fig. 3. Results of the position measurements of the maser feature in Orion KL. (a) The movement of the mascr fea- 
ture in right ascension as a function of time. (b) The same as (a) in dechnation. (c) The movement of the 
maser feature on the sky. Solid lines represent the best fit model with the annual parallax and linear proper mo- 
tion for the maser feature. Filled circles represent the observed positions of the maser feature with error bars indicat- 
ing the standard deviations of the least-squares analysis as listed in Table 1 (0.36 mas in right ascension and 0.74 mas 
in declination). The reference position is the same as in Figure 2. Observed epochs are indicated in the panel (c). 



than due to an inappropriate model in the least-squares 
analysis. In this section, wc will consider possible sources 
of these astrometric errors. 

As reported previously in the literature (Kurayama et 
al. 2005; Hachisuka et al. 2006; Honma et al. 2007; 
Sato et al. 2007), it is difficult to estimate the individual 
error sources in the VLBI astrometry quantitatively. We 
therefore estimate the uncertainties in the measured posi- 
tion of the maser feature to be 0.36 mas and 0.74 mas in 
right ascension and declination, respectively, based on the 
standard deviations of the least-squares analysis as listed 
in Table 1. The standard deviations obtained in this pa- 
per are larger than those of previous observations with 
VERA (Honma et al. 2007; Sato et al. 2007), especially 
in declination. 

The most serious error source in the VLBI astrome- 
try in the 22 GHz band is likely to be the atmospheric 
zenith delay residual due to the troposphcric water va- 
por. This is caused by the difference in the optical path 



lengths through the atmosphere between the target and 
reference sources because the elevation angle of the tar- 
get source is usually different from that of the reference 
source. According to the discussions in Honma et al. 
(2007), a path length error due to the atmospheric zenith 
delay residual of 3 cm would cause a position error of 0.04- 
0.12 mas in the case of a separation angle between the tar- 
get and reference sources of 0.7 degrees at the elevation 
angle of 20-90 degrees. If we consider an extreme example, 
with the observed elevation angle of 20 degrees and the 
atmospheric zenith delay residual of 10 cm, the position 
error in the observations of Orion KL and J0541— 0541, 
with a separation angle of 1.62 degrees, is estimated to be 
0.75 mas. This value is clearly overestimated because the 
path length errors should be suppressed at the higher ele- 
vation angle. Furthermore, the atmospheric zenith delay 
residual of 10 cm is unrealistic because we have corrected 
such a large residual before phase-referencing imaging. 
Therefore, the atmospheric zenith delay residual alone 
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Table 2. Results of the proper motion measurements for the observed maser feature and source I 







Absolute proper motion 




Proper 


motion relative to source I 


Source 

Name 


Ma 

(mas yr""'^) 


M.5 A* 
(mas yr""'^) (mas yr~^) 


Vt 

(km s~^) 


IJ-l 4 
(mas yr""'^) (mas yr""'^) (mas yr^"'^) (km s""'^) 


Maser" 
source P 


2.77(0.09) 
3.5(0.7) 


-8.97(0.21) 9.39(0.20) 
-4.4(0.7) 5.6(0.7) 


19.7(0.4)" 
12(2)'' 


-0.7(0.7) 
0.00 


-4.6(0.7) 4.6(0.7) 9.7(1.5)" 
0.00 0.00 0.00 



Note — Numbers in parenthesis represent the estimated uncertainties. 

a: Absolute proper motion of source I is taken from Rodriguez et al. (2005). 

h: Calculated assuming the distance of 437 pc. 



cannot fully explain our position errors, although it would 
contribute to the large part of the error source in our as- 
trometry, especially in declination. 

One of other possibilities for the error sources in the 
observed position is a variation of the structure in the 
maser feature. With regard to this, we confirmed that 
peak positions of the maser spots within the analyzed 
maser feature were sometimes shifted by about 0.2 mas 
from those of the adjacent channels. In addition, the sys- 
tematic velocity shift from 25.7 km s~^ to 24.5 km s~^ 
was observed during the observing period of 2 years, indi- 
cating the variation of the maser feature. Although there 
is no reason that the structure in the maser feature affects 
the astrometric accuracy only in declination, it would be 
one of the major sources of errors in the astrometry with 
the H2O maser lines as well as the atmospheric zenith 
delay residual. The effect of the spatial structure of the 
maser feature is more significant for Orion KL than the 
other sources (Kurayama et al. 2005; Hachisuka et al. 
2006; Honma et al. 2007; Sato et al. 2007) because the 
distance to Orion KL (437 pc) is nearer than the others 
by a factor of 5-10 (2-5 kpc). However, this effect is in- 
versely proportional to the distance to the target source 
just the same as its annual parallax. This means that 
the annual parallaxes of the more distant sources can be 
measured with almost the same precision as in the case 
of Orion KL, if the dominant error source in astrometry 
is due to the structure effect rather than the atmospheric 
zenith delay residual. In fact, the relative uncertainty in 
the annual parallax of the further source, S269, is found 
to be comparable to that of Orion KL, about 4%, in the 
case of using the data for right ascension only (Honma et 
al. 2007). Further VLBI observations of maser features 
with shorter baselines should be able to confirm this effect, 
with which more extended structures of maser features are 
imaged. 

On the other hand, the variation of the structure of 
the reference source, J0541— 0541, would be negligible for 
the measurements of the annual parallax and proper mo- 
tion because we found no evidence for significant structure 
of J0541— 0541 in our observations. The uncertainty in 
the absolute position of the reference source J0541— 0541, 
0.28 mas and 0.46 mas in right ascension and declination, 
respectively (Ma et al. 1998), also does not affect the 



derived annual parallax and proper motion because this 
uncertainty gives only a constant offset to the position of 
the maser feature. According to the discussions in Honma 
et al. (2007), astrometric errors in the VERA observa- 
tions arising from uncertainties in the station position, 
delay model, and path length errors due to ionosphere are 
estimated to be smaller by an order of magnitude, and 
hence, they do not have significant effects on astrometric 
accuracy. Therefore, we conclude that the major sources 
of our astrometric errors are due to the atmospheric zenith 
delay residual and variability of the structure of the maser 
feature. 

4-.2. Annual Parallax and Distance to Orion KL 

We successfully obtained the annual parallax of 
Orion KL to be 2.29±0.10 mas, corresponding to the dis- 
tance of 437±19 pc. This is the first time that the distance 
to Orion KL is determined based on the annual parallax 
measurements. Genzel et al. (1981) derived the dis- 
tance to Orion KL to be 480 ± 80 pc from the statistical 
parallax method, using proper motions and radial veloc- 
ities of the H2O maser features. Our result is consistent 
with that of Genzel et al. (1981), although the accu- 
racy of our measurements is significantly improved. The 
most important progress in our new results is due to the 
geometric nature of our measurements without any as- 
sumption unlike the statistical parallax method, in which 
appropriate kinematic modeling for Orion KL is required 
(Genzel et al. 1981). The accuracy of the annual par- 
allax measurements in our study is limited mainly due 
to the atmospheric zenith delay residual and the struc- 
ture of the maser feature, both of which are difficult to be 
predicted and measured completely in the current obser- 
vational study. In principle, it will be possible to achieve 
much higher precision using the results of all the maser 
features in Orion KL, which will reduce the statistical er- 
ror by a factor of N~^'^ where N represents the number 
of observed maser features. This expectation will be con- 
firmed in the further analyses of the VERA observations. 

4.3. Absolute Proper Motion of the Maser feature in 
Orion KL 

Along with the annual parallax measurements, we 
successfully detected the absolute proper motion in 
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our phase-referencing astrometry with VERA. Figure 2, 
Tables 1 and 2 show the absolute proper motion of the 
maser feature in Orion KL. At the distance of 437 pc, the 
proper motion of 1 mas yr"-'^ corresponds to the transverse 
velocity of 2.1 km s^^. The observed absolute proper 
motion of the H2O maser feature (2.77±0.09 mas yr~^ 
and — 8.97±0.21 mas yr~^ in right ascension and declina- 
tion, respectively) corresponds to 9.39±0.20 mas yr^^ or 
19.7±0.4 km s"^ toward south. 

Recently, Rodriguez et al. (2005) and Gomez et al. 
(2005) measured the proper motion of radio continuum 
sources in the Orion KL region with the VLA, as shown 
in Figure 2 and Table 2. Subtracting the proper motion 
vector of source I from that of the observed maser feature, 
we can obtain the proper motion of the maser feature with 
respect to source L As Gomez et al. (2005) have already 
mentioned, the precision of the absolute proper motion 
measurements by Rodriguez et al. (2005) is higher than 
that by Gomez et al. (2005). Therefore, we adopt the 
proper motion of source I inferred by Rodriguez et al. 
(2005), 3.5±0.7 mas yr~^ and — 4.4±0.7 mas yr~^ in right 
ascension and declination, respectively, in the following 
discussions. The proper motion of the maser feature with 
respect to source I is inferred to be — 0.7±0.7 mas yr~^ 
and — 4.6±0.7 mas yr"-'^ in right ascension and declination, 
respectively, as listed in Table 2. The magnitude of the 
proper motion is 4.6±0.7 mas yr~^ or 9.7±1.5 km s~^ 
toward south with a position angle of —171 degrees, which 
agrees well with the direction of the outflow powered by 
source I. Therefore, we conclude that the absolute proper 
motion of the observed maser feature is the sum of outflow 
motion powered by source I and the systematic motion of 
source I itself. 

However, a detailed model of the outflow powered by 
source I is still debatable. Greenhill et al. (1998) 
first proposed that the biconical high-velocity outflow 
traced by the SiO maser lines lies along the northwest- 
southeast direction, while the low-velocity equatorial out- 
flow traced by the H2O maser lines exists along the 
northeast-southwest direction. On the other hand, they 
changed the interpretation based on the recent results 
that the outflow is along the northeast-southwest direc- 
tion, which is perpendicular to the first model, and that 
the SiO maser lines trace the edge-on disk perpendicu- 
lar to the outflow (Greenhill et al. 2004). We cannot 
distinguish these two different models in this paper be- 
cause the distribution of the H2O masers, elongated along 
the northeast-southwest direction as shown in Figure 2, is 
consistent with both models and in addition, the proper 
motion of the observed H2O maser feature is almost in- 
termediate (toward south) between the proposed outflow 
axes (Greenhill et al. 1998, 2004). The velocity structure 
in the Orion KL region is quite complicated as Greenhill 
et al. (2004) suggested, and hence, further discussions 
about the proper motions of all the H2O maser features 
arc required to construct the detailed model of the out- 
flow in the Orion KL region, which will be presented in a 
forthcoming paper. 
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